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Phytoplankton productivity across prairie saline 
lakes of the Great Plains (USA): a step toward 
deciphering patterns through lake classification 
models 

Courtney R. Salm, Jasmine E. Saros, Sherilyn C. Fritz, Christopher L. Osburn, and 
David M. Reineke 


Abstract: We investigated patterns of primary production across prairie saline lakes in the central and northern Great 
Plains of the United States. Based on comparative lake sampling in 2004, seasonal predictors of algal primary productivity 
were identified within subsets of similar lakes using a combination of Akaike’s information criterion (AIC) and classifica¬ 
tion and regression trees (CART). These models indicated complex patterns of nutrient limitation by nitrogen (N), phos¬ 
phorus (P), and iron (Fe) within different lake groups. Nutrient enrichment assays (control, + Fe, + N, + P. + N + P) were 
performed in spring and summer of 2006 to determine if phytoplankton in selected lakes followed predicted patterns of nu¬ 
trient limitation. Both the comparative lake sampling and experimental results indicated that N limitation was widespread 
in these prairie lakes, with evidence for secondary P limitation in certain lakes. In the experiments, iron did not stimulate 
primary production. Our results suggest that given the diverse geochemical nature of these lakes, classification models that 
separate saline lakes into subsets may be an effective method for improving predictions of algal production. 

Resume : Nous avons etudie les patrons de production primaire dans des lacs sales de prairie repartis dans le centre et le 
nord des Grandes Plaines des Etats-Unis. A partir d’un echantillonnage comparatif des lacs en 2004, nous avons identifie 
les variables predictives saisonnieres de la production primaire des algues dans des sous-ensembles de lacs en utilisant 
conjointement le critere d’information d’Akaike (AIC) et des arbres de classification et de regression (CART). Ces mod¬ 
ules identifient des patrons complexes de limitation par les nutriments d’azote (N), de phosphore (P) et de fer (Fe) dans 
les differents groupes de lacs. Nous avons fait des tests d’enrichissement de nutriments (temoin, + Fe, + N, + P, + N + P) 
au printemps et a Fete 2006 afin de determiner si le phytoplancton dans les lacs selectionnes suivait les patrons predits de 
limitation par les nutriments. Tant Fechantillonnage comparatif des lacs que les resultats experimentaux indiquent que la 
limitation par N est commune dans ces lacs de prairie, avec des indications d’une limitation secondaire par P dans certains 
lacs. Dans les experiences, le fer ne stimule pas la production primaire. Nos resultats indiquent qu’etant donne la nature 
geochimique diverse de ces lacs, les modeles de classification qui separent les lacs sales en sous-ensembles peuvent etre 
des methodes efficaces pour ameliorer les predictions de production des algues. 

[Traduit par la Redaction] 


Introduction 

Prairie saline lakes are widely distributed across semi-arid 
regions of the world (Williams 1981) and are important hab¬ 
itats for migratory waterfowl and carbon sequestration (Batt 
et al. 1989; Euliss et al. 2006). Fossil algal records from 
these lakes are often used to reconstruct drought frequency 


and severity over time scales ranging from centuries to mil¬ 
lennia (Cumming and Smol 1993; Laird et al. 1998; Fritz et 
al. 2000). However, although some of these lakes are limited 
by phosphorus or trace metals (Waiser and Robarts 1995; 
Evans and Prepas 1997), broad trends in the factors that 
control algal production across suites of these systems are 
not yet apparent (Bierhuizen and Prepas 1985; Campbell 
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and Prepas 1986). Commonly, these lakes do not follow the 
paradigm established for freshwater lakes (Waiser and Ro- 
barts 1995) in which spring total phosphorus (TP) concentra¬ 
tions predict summer chlorophyll a (Chi a) concentrations 
(Dillon and Rigler 1974; Smith 1979). In prairie saline 
lakes, TP concentrations are often greater than 50 pg-L ', 
and total nitrogen (TN) concentrations can exceed 
1000 pg-L -1 (Campbell and Prepas 1986). Despite these 
high concentrations of nutrients, the algal biomass is lower 
than predicted based on freshwater models (Campbell and 
Prepas 1986; Robarts et al. 1992; Evans and Prepas 1997). 

Several factors may account for the poor performance of 
the freshwater spring TP - summer Chi a model in saline 
lakes. Although total nutrient concentrations are high, a sub¬ 
stantial fraction of the nutrients may not be bioavailable 
(Waiser and Robarts 1995), particularly because concentra¬ 
tions of dissolved organic material (DOM) are very high in 
these lakes, ranging from 20 to 800 mg C L 1 (measured as 
dissolved organic carbon (DOC); Curtis and Adams 1995; 
Evans and Prepas 1997; Arts et al. 2000). DOM can com¬ 
plex with nutrients and alter their bioavailability (Hessen 
and Tranvik 1998; Bushaw-Newton and Moran 1999; Fin¬ 
dlay and Sinsabaugh 2003). In addition, the ionic composi¬ 
tion and high salinity of saline lakes affects nutrient cycling 
and availability, as well as uptake by algae (Caraco et al. 
1989; Waiser and Robarts 1995; Saros and Fritz 2000). 
Chlorophyll concentrations may be affected by intense zoo¬ 
plankton grazing on standing phytoplankton crops in saline 
lakes (Anderson 1958), and limited information is available 
on rates of primary production rather than chlorophyll con¬ 
centrations alone (Armstrong et al. 1966; Hammer 1981). 
These conditions suggest that measuring total and dissolved 
nutrient pools may not provide enough information to dis¬ 
cern the factors that drive algal production in these lakes; 
physiological parameters of algal production and nutrient 
status may be more useful in clarifying any patterns (Waiser 
and Robarts 1995). 

An additional factor that may hinder the development of a 
predictive model of phytoplankton production in saline lakes 
is the assumption that one model and one variable can have 
sufficient explanatory power to model phytoplankton re¬ 
sponse in these chemically complex systems. Unlike marine 
systems, the dominant ions in saline lakes vary extensively 
(Fritz et al. 1993; Gosselin 1997), and the cycling of phos¬ 
phorus and nitrogen can differ greatly, for example, between 
bicarbonate- and sulfate-dominated systems (Cole et al. 
1986; Caraco et al. 1989). These differences suggest the 
need to develop a classification system to separate these di¬ 
verse lakes into subsets based on similar characteristics. 
Lake classification has been used for a variety of applica¬ 
tions, from assessing trophic state to managing water quality 
(Emmons et al. 1999; Spndergaard et al. 2005; Bulley et al. 
2007). 

We examined patterns of primary production in phyto¬ 
plankton across a suite of lakes in the central and northern 
Great Plains of the United States in the spring and summer 
of 2004. In addition to measuring environmental parameters, 
physiological methods were used to assess the nutrient status 
of phytoplankton by measuring alkaline phosphatase activity 
and seston nutrient ratios. Rather than employing chloro¬ 


phyll measurements alone as in other studies, we directly 
measured phytoplankton primary production across all lakes 
via 14 C uptake rates. Subsequently, we tested the model pre¬ 
dictions from the survey data in Cubitainer experiments in a 
set of lakes in the spring and summer of 2006. 

Materials and methods 

Study sites 

Prairie saline (>3 g-L *) and subsaline (0.5-3 g-L 1 ) lakes 
included in this study are located in the central (Nebraska) 
and northern (North Dakota, South Dakota, and Montana) 
Great Plains (CGP and NGP, respectively) (Fig. 1). On aver¬ 
age, the Great Plains are characterized as semi-arid with a 
negative effective moisture balance (precipitation minus 
evapotranspiration, P - E). Land-use patterns vary across 
the region and include primarily cropland, rangeland, and 
undisturbed native grassland. 

In the NGP, thousands of lakes occur east of the Missouri 
River; most of these are morainal depressions or dammed 
river valleys formed upon retreat of the Pleistocene ice mar¬ 
gin (Bluemle and Clayton 1984). In eastern and central 
North Dakota and South Dakota, the majority of saline lakes 
are sulfate-dominated as a result of the widespread presence 
of pyrite in bedrock and surface deposits (Fritz et al. 1993). 
Major cation concentrations and proportions are highly vari¬ 
able. A smaller number of lakes are found in the glaciated 
regions of northeastern Montana - northwestern North Da¬ 
kota immediately north of the Missouri River. The majority 
is in paleochannels of the Missouri River system and is do¬ 
minated by carbonates because of inflow from shallow 
groundwater systems in calcareous outwash and till (Do¬ 
novan 1994). 

Lakes are also widespread in the Nebraska Sandhills of 
the CGP, a large region of Holocene eolian sands overlying 
Pleistocene and late-Tertiary alluvial sands and silts, which 
is the principal recharge area of the High Plains (Ogallala) 
Aquifer. Here lakes are formed by deflation and blockage 
of interdune valleys by eolian sands, and most lakes are 
thought to have formed during dry intervals during the Hol¬ 
ocene (Loope et al. 1995; Mason et al. 1997). The majority 
of lakes is dominated by bicarbonate, although a small num¬ 
ber of sulfate-enriched lakes occur. Sodium and potassium 
are the most abundant cations (LaBaugh 1986; Gosselin 
1997). 

Lakes across the NGP and CGP are topographically 
closed, but many are hydrologically connected to ground- 
water. Salinities range from 0.1 to greater than 100 g-L 1 
(note that for simplicity, we use the term “prairie saline 
lakes” to refer to the entire group of lakes, recognizing that 
technically some of the lakes are subsaline). Lakes were se¬ 
lected to maximize variation in conductivity, ion composi¬ 
tion, and nutrient concentrations. In the spring, these lakes 
were primarily dominated by diatoms (largely Fragilaria 
crotonensis and Cyclotella quillensis/meneghiniana) with 
some cyanobacteria ( Gloeocapsa sp. and Aphanocapsa sp.). 
In the summer, cyanobacteria dominated the phytoplankton 
assemblages (Aphanizomenon sp. and Aphanocapsa sp.) 
with some diatoms ( Cyclotella quillensis and Surirella sp.) 
(Salm et al. 2009). 
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Fig. 1 . Map of central Great Plains, USA. Ellipses indicate study 
sites for the 2004 survey. NGP, northern Great Plains; CGP, central 
Great Plains. 
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Lake sampling and model development 
Data collection 

Data for approximately 30 parameters were collected for 
lakes across the study area in both spring (24 lakes) and 
summer (30 lakes), including temperature, conductivity, pH, 
total alkalinity, ion composition, Chi a, nutrients (total, par¬ 
ticulate, and dissolved), DOM characteristics, alkaline phos¬ 
phatase activity (APA), and rates of primary production. For 
these analyses, approximately 10 L of lake water were col¬ 
lected with a van Dorn horizontal bottle from a depth of 1 m 
in lakes that were at least 2 m deep and from the midpoint 
of the water column in lakes that were <1 m deep. 

Temperature and conductivity were measured with a port¬ 
able conductivity meter (WTW MultiLine P4); pH was meas¬ 
ured with a pH meter (Corning); and total alkalinity was 
determined by titration (American Public Health Association 
(APHA) 1998). Cations were measured via atomic absorption 
spectroscopy (Varian 220FS with a GTA-110 graphite fur¬ 
nace and a VGA-77 vapor generation unit), and anions were 
measured by ion chromatography (Dionex ICS-90). 

In the field, whole water samples were collected for total 
P and total N and acidified with H 2 SO 4 . Samples for total 
dissolved P, soluble reactive P, nitrate + nitrite, and dis¬ 
solved Si were filtered through 0.45 |im Millipore mem¬ 
brane filters. Samples for total P, dissolved P, and 
particulate P were first digested with potassium persulfate 
and measured by ascorbic acid methods (Lind 1985; APHA 
1998), as were those for soluble reactive P. Total N was 
measured by alkaline potassium persulfate digestion (D’Elia 
et al. 1977) and the ultraviolet (UV) absorption method 
(APHA 1998), and nitrate + nitrite N was determined by 
the hydrazine reduction method (Downes 1978). Dissolved 
Si was measured following the methods of Wetzel and Lik¬ 
ens (1991). Water was also filtered through 0.7 pm glass- 
fiber filters (Whatmann GF/F) for analysis of Chi a, as 
well as particulate carbon and nitrogen analyses (in dupli¬ 
cate). Filters were collected onto Petri dishes and either 


wrapped in foil and frozen (Chi a) or refrigerated (particu¬ 
late C and N) until processing. Chlorophyll was analyzed 
spectrophotometrically after pigment extraction with 90% 
acetone (Varian Cary-50 UV-VIS Spectrophotometer; 
APHA 1998) within two weeks of collection to maintain 
sample integrity. Filters for particulate C and N were 
fumed with concentrated HC1 in a glass desiccator and 
measured by combustion and gas chromatography with an 
elemental analyzer (Costech). 

Dissolved iron was quantified in samples that were fil¬ 
tered through 0.45 pm Millipore membrane filters and acidi¬ 
fied at the time of collection. Total recoverable iron was 
determined on separate acidified samples prepared by mix¬ 
ing unfiltered lake water with hydrochloric acid and nitric 
acid to improve the solubility of iron during digestion. Sam¬ 
ples were covered and digested overnight (~12 h) at 90 °C 
on a hotplate and diluted in 1 % nitric acid prior to analysis. 
Samples for total and dissolved iron were analyzed with an 
inductively coupled plasma mass spectrometer (GV Instru¬ 
ments Platform XS). Total iron was quantified based on the 
response of 56 Fe. Potential interferences were monitored in 
samples and blanks fortified with these elements, and In¬ 
dium was added to all samples and standards to monitor for 
additional matrix effects. The method detection limit was 
1.0 ppb. 

Samples for DOM characterization were prefiltered 
through 0.2 pm pore-size membranes before analysis. DOC 
concentration (mg C L ' ) was measured by wet chemical ox¬ 
idation on an OI Analytical 1010 TOC analyzer, following 
the recommendations of Osburn and St. Jean (2007) for 
high-salinity samples. Sodium persulfate (450 g-L _1 ; cleaned 
by heating to a near-boil and then rapidly cooling) was 
added to the reactor and allowed to react for 10 min, con¬ 
verting all DOC to CO 2 . The CO 2 was quantified by nondis- 
persive IR detection and calibrated to potassium biphthalate 
standards over the range of 1-100 mg CL '. Sample vol¬ 
umes of 100-2000 |xL were injected to stay within this cali¬ 
bration range. 

Absorption of water samples was measured from 250 to 
650 nm on a dual-beam spectrophotometer (Shimadzu 
UV1601) against Milli-Q water in 1 cm cuvets due to the 
strong absorbance of these samples. Raw absorbance was 
converted to absorption coefficients: 

a(X) = A(X) x 2.303/D 

where A(X) is the raw absorbance of the sample at wave¬ 
length X, and D is the pathlength of the 1 cm cuvet (in 
metres) (Kirk 1994). The constant of 2.303 converts from 
the natural logarithm. The ratio of absorption at 250- 
365 nm has been used to observe relative changes in mole¬ 
cular weight of DOC (De Haan and De Boer 1987). An in¬ 
crease in the ratio indicates a decrease in molecular weight. 

Bulk alkaline phosphatase activity was measured in the 
field with the 4-methylumbelliferyl phosphate method (Hill 
et al. 1968) on a field fluorometer (Turner Designs Model 
10-AU Field Fluorometer). Bulk APA:Chl a ratios were 
used to normalize enzyme activity to the amount of algal 
biomass. 

Rates of primary productivity in each lake were assessed 
using a modified light-dark bottle 14 C uptake assay (Wetzel 
and Likens 1991). Bottles (300 mL) were incubated under a 
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suspended light bank (Freshwater AquaLight by Coralife, 
-450 p/E-cm 2 -s r ) in a water bath held at a temperature con¬ 
sistent with ambient seasonal lake conditions (2004: spring, 
15 °C; summer, 20 °C). We chose this approach, rather than 
in situ incubations in each lake, to maintain similar light and 
temperature conditions across all lakes. Bottles were filled 
with lake water prescreened through 212 |xm Nitex mesh to 
remove large grazers. During the summer, when colonial cy¬ 
anobacteria populations were present in lakes, prefiltering 
was not possible for some lakes. In these instances, grazers 
were handpicked from bottles with a pipette. Each bottle 
was inoculated with 0.5-2.0 mL NaH 14 C 03 solution (Perkin 
Elmer, 1 trCi-mL 1 ) and incubated onshore for 3 h at some 
point between 0900 and 1500 h (within the diurnal period of 
maximum algal productivity). After incubation, two aliquots 
(the size of which depended on algal densities) were col¬ 
lected from each bottle on 0.45 |im Millipore HA filters, 
which were dried in a desiccator overnight and stored in a 
cool, dry location until processing in the laboratory. Filters 
were fumed with concentrated HC1 in a glass desiccator to 
remove residual inorganic 14 C and carbonates. They were 
dissolved in 10 mL of Filter-Count scintillation cocktail 
(Perkin Elmer) and read on a scintillation counter (Beckman 
Coulter LS 6500; 32 768 channels with 0.06 keV per chan¬ 
nel resolution), with primary production rates calculated ac¬ 
cording to Wetzel and Likens (1991). 

Data analysis 

Values for all variables except pH, seston ratios, and 
a 250 :ci 365 were log 10 -transformed to correct for unequal var¬ 
iance after Levene’s homogeneity of variance test, and sim¬ 
ple and stepwise multiple regression were used to identify 
predictors of primary production (a = 0.05). In addition to 
these methods, Akaike’s information criterion (AIC; Burn¬ 
ham and Anderson 2002) was used to explore how multiple 
factors may control primary production across all lakes in 
the data set. Although a priori hypotheses are strongly sug¬ 
gested when using AIC, in cases such as this when it is dif¬ 
ficult to formulate reasonable a priori models, AIC may be 
used in conjunction with secondary studies to confirm im¬ 
portant predictor variables (Burnham and Anderson 2002). 
A program for running AIC analyses in SAS (version 9.1; 
SAS Institute Inc., Cary, North Carolina) was written by 
Episystems, Inc. (St. Paul, Minnesota) and used to analyze 
the spring and summer comparative lake sampling data. The 
best set of approximating models was found from a larger 
set of candidate models by selecting those with the lowest 
AIC. It has been suggested that models for which the AIC 
value is within 4 units of the lowest would constitute an ap¬ 
proximate 95% confidence set of top models (Burnham and 
Anderson 2002); thus this criterion was used to determine 
which models were retained, parameter weights, and model- 
averaged parameter estimates. 

Exploratory classification models also were developed 
with classification and regression trees (CARTs; De’ath and 
Fabricius 2000) using the program R (version 2.1.1; http:// 
www.r-project.org/), with the rate of primary production as 
the predictor variable for regression trees and phosphorus 
limitation status (severely limited or not limited) as the pre¬ 
dictor variable for classification trees. An APA : Chi a value 


above 0.005 is considered indicative of severe P limitation 
(Healey and Hendzel 1980); hence this was used as the di¬ 
viding value between lakes that were potentially P-limited 
(APA : Chi a > 0.005) or not (APA : Chi a < 0.005). Simple 
regressions were then performed in SPSS (version 11.5 for 
Windows; SPSS Inc., Chicago, Illinois) to identify parame¬ 
ters controlling primary production within each of the result¬ 
ing groups. 

Experiments and model testing 
Experimental design 

Using the groups identified in the CART analyses (see 
Results below), nutrient enrichment experiments were con¬ 
ducted in spring and summer 2006 as a test of these models. 
Assays were similar in spring and summer, with five treat¬ 
ments created with the following additions: control (no nu¬ 
trient addition), Fe, N, P, and N + P (n = 3). For the spring 
and summer experiments, a subset of two or three lakes was 
chosen from within each group identified by the CART 
analyses. Water from each lake was filtered through a 
212 p,m mesh to remove large zooplankton and incubated 
in 4 L Cubitainers (VWR TraceClean™ Cubitainer™ con¬ 
tainers, made of low-density polyethylene) in Coldwater 
Lake at approximately 0.5 m depth in the spring (15 °C) 
and just below the surface in the summer (20 °C). Appropri¬ 
ate iron (11.7 pmol Fe L 1 in the form of FeCl3-6H20, added 
along with 11.7 pmol EDTAL 1 ), nitrogen (18 pmol N-L ! 
in the form of NaNCy, and phosphorus (5 pmol P L -1 in 
the form of NaHiPCF) additions were made to each Cubi¬ 
tainer. Transmission scans of Cubitainer plastic on a spectro¬ 
photometer (Varian Cary-50 UV-VIS Spectrophotometer) 
indicated that for the visible spectrum, approximately 60%- 
75% of ambient light was allowed through the plastic, 
whereas at lower wavelengths (-300 nm), 40%-50% of am¬ 
bient light was allowed through. For each Cubitainer, rates 
of primary production ( 14 C uptake assay) were measured in¬ 
itially and 3 days after nutrient additions, because prelimi¬ 
nary laboratory experiments indicated that, based on daily 
measurements over a 7-day period, the greatest changes in 
production rates occurred on day 3 (data not shown). One 
light bottle was analyzed for each Cubitainer, and two dark- 
bottle replicates were analyzed for each lake. The dark-bottle 
replicates were taken from a random control Cubitainer and a 
random N + P Cubitainer for each lake to assess whether the 
higher biomass in nutrient addition treatments affected the 
14 C measurements. Previous use of this method during the 
comparative lake sampling in 2004 showed very little varia¬ 
tion among dark-bottle replicates. 

Data analysis 

Results from these experiments were statistically analyzed 
using analysis of variance (ANOVA) in SPSS (version 11.5 
for Windows; SPSS Inc., Chicago, Illinois) to determine if 
rates of primary production were significantly different 
among treatments (a = 0.05). Levene’s homogeneity of var¬ 
iance test was used to check for equal variances, and rates 
of primary production were logio-transformed to correct for 
unequal variances. Tukey’s post-hoc analysis was also used 
to compare mean values across treatments. 
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Results 

Lake sampling and model development 

Rates of primary production were generally higher across 
lakes in summer compared with spring (Table 1). Regional 
differences among lakes also were evident, as CGP lakes 
had much higher rates of productivity than those of the 
NGP based on the 14 C uptake assays (spring average rates 
of primary productivity: CGP, 263.4 mg C-m 3 * NGP, 
33.0 mg C-m 3 * *; summer average rates: CGP, 526.5 mg 
C-m- 3 *- 1 ; NGP, 124.9 mg C-m- 3 *- 1 ). CGP lakes had high 
concentrations of nutrients and DOC and were shallow 
(<1 m) and very turbid, making filtration difficult (only pos¬ 
sible to filter 1-2 mL at most through 0.45 pm filters). 
Some of these systems dried out in the summer, preventing 
parallel seasonal studies. Given these difficulties with the 
CGP lakes, we examined patterns in the full data set of 
lakes, as well as the subset of lakes (n =18) found in the 
NGP alone. 

Simple regressions between rates of primary production 
and key nutrients did not indicate strong limitation by a sin¬ 
gle nutrient (Table 2), with the exception of N 03 “ for all 
lakes in the spring ( R 2 = 0.519, F = 28.08, df = 1,22, p < 
0.001). Other low R 2 but significant simple regressions with 
nutrient parameters included total P (R 2 = 0.306, F = 9.718, 
df = 1,22, p = 0.005) and total N (R 2 = 0.202, F= 5.561, df = 
1,22, p < 0.028) for all lakes in the spring, total P (R 2 = 
0.316, F = 12.951, df = 1,28, p < 0.001) and total N (R 2 = 
0.168, F = 5.634, df = 1,28, p = 0.025) for all lakes in the 
summer, and total N (R 2 = 0.238, F = 4.985, df = 1,16, p = 
0.040) for NGP lakes only. 

Multiple regressions revealed variable predictors of pri¬ 
mary production by season and region. For the spring, NO 3 - 
was the only variable selected by stepwise regression for all 
lakes (R 2 = 0.519, F = 22.66, df = 1,21, p < 0.001), with a 
positive relationship between this parameter and primary 
production rates. This relationship was not found for the 
subset of NGP lakes alone, in which no variables were sig¬ 
nificant predictors of primary production. In the summer, 
stepwise selection indicated C:P and dissolved Fe as the 
best predictors for all lakes ( R 2 = 0.606, F = 20.75, df = 
1,28, p < 0.001), and absorptivity at 350 nm (indicative of 
chromophoric dissolved organic matter) for NGP lakes alone 
(R 2 = 0.311, F = 7.23, df = 1,16, p = 0.016). These parame¬ 
ters were positively correlated to rates of primary produc¬ 
tion. 

AIC models were developed separately for spring and 
summer; within each season, models for all lakes, as well 
as just the NGP lakes, were generated (Table 3). All of the 
models suggested multivariate control over primary produc¬ 
tion across these lakes. Each model also identified one or 
two of the strongest variables (with a predictor weight close 
to or equal to 1 ) correlated with rates of primary production. 
Looking across the variables with the highest predictor 
weights, the AIC approach identified dissolved Fe and NO 3 - 
as the best predictors of primary production in spring for all 
lakes, and calcium, dissolved Fe, and TP for the NGP lakes 
subset. In the summer, the best predictors for all lakes were 
C:P and ions (CF, K + , Na + ), and phosphorus parameters 
(TP, N:P, and SRP) for the NGP lakes. Thus, results indi¬ 


cated that the limiting factors for rates of primary produc¬ 
tion differed regionally as well as seasonally. 

In the CART analyses, spring regression trees showed a 
split, with NO 3 - as the only branching variable (Fig. 2). This 
separation appeared to be related to geographic distribution, 
as high nitrate lakes ([NO 3 -] > 17.5 pg-L ' ) were located in 
the CGP and low nitrate lakes ([N 03 “] < 17.5 pg-L 1 ) were 
located in the NGP. Classification further split the NGP 
lakes by APA: Chi a ratio. Within the three groups de¬ 
lineated in Fig. 2, simple regressions did not yield any sig¬ 
nificant predictor variables for rates of primary production 
at the a = 0.05 level. However, regressions for high NO 3 - 
lakes from the CGP showed a positive correlation between 
primary production and maximum depth (R 2 = 0.5000, F 
Till. df = 1,4, p = 0.12), which may be a reflection of the 
very high turbidity and greater potential light limitation in 
the shallower lakes of the CGP. Low NO 3 -, low APA: Chi a 
lakes (designated SPR-N to reflect potential N limitation) 
showed a positive correlation between primary production 
and N:P ratios (R 2 = 0.3623, F = 3.41, df = 1,6, p = 0.11), 
and low NO 3 -, high APA: Chi a lakes (designated SPR-NP 
to reflect potential N and P co-limitation) showed a positive 
correlation between primary production and unfiltered APA 
rates (R 2 = 0.3323, F = 3.98, df = 1,8, p = 0.08) such that 
lakes in this subset with higher APA rates were more pro¬ 
ductive. 

Summer regression trees showed a split with C:P seston 
ratios (Pig. 3). The high C:P group (C:P > 3591) contained 
most of the highly productive CGP lakes, whereas the group 
with lower C:P values (C:P < 3591) had lower rates of pri¬ 
mary production on average and included NGP lakes and the 
remainder of the CGP lakes. The split was at an extremely 
high C:P ratio, outside of typical ranges for seston ratios for 
freshwater lakes (Elser et al. 2000; Sterner and Elser 2002), 
and generally follows geographic distribution as in the 
spring. When CGP lakes were excluded from the regression 
tree, there were no branches within the NGP lakes. These 
data suggested that phosphorus gradients were involved in 
controlling primary production during the summer, so 
APA: Chi a classification trees were used again for the 
NGP lakes alone. The main branching factor in this subset 
of lakes was soluble reactive P (SRP), with a branch point 
of 25 pg-L 1 . Simple regressions for the higher SRP group 
(designated SUM-N) showed a positive correlation between 
primary production rates and N:P ratios (R 2 = 0.38, F = 
5.55, df = 1,9, p = 0.043). The regression for the lower SRP 
group (designated SUM-P) between rate of primary produc¬ 
tion and total P was significant (R 2 = 0.78, F = 17.5, df = 
1,5, p = 0.009). 

Spring experimental results 

Lor the spring experiment, three lakes were selected from 
the SPR-NP group (Alkaline Lake, Coldwater Lake, and 
Lree People Lake), as well as two from the SPR-N group 
(East Devils Lake and Stink Lake). The nutrient treatments 
affected primary production rates in four of the five lakes 
tested (Lig. 4; Table 4). Of the SPR-NP lakes predicted to 
respond to both N and P, rates in Alkaline Lake (F = 30.6, 
df = 4,10, p < 0.001) and Coldwater Lake (F = 75.6, df = 
4,10, p < 0.001) increased in the N treatment, with addi¬ 
tional increases in the N + P treatments. Lree People Lake 
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Table 1 . Summary of 2004 spring and summer survey data used in analyses. 


Parameter 

Units 

All lakes 

CGP only 

NGP only 

Spring 

Primary production rate 

mg C-m a -h 1 

91 (4—452) 

263 (138-452) 

33 (4-102) 

Chl a 

Pg'L' 

76 (1-1060) 

225 (2-1060) 

10 ( 1 - 22 ) 

7 

max 

m 

4.1 (1.0-13.0) 

1.6 (1-2.5) 

5.1 (1.0-13.0) 

Temperature 

°c 

15 (9-22) 

18 (16-22) 

13 (9-17) 

pH 


9.04 (7.94-10.43) 

9.61 (8.89-10.43) 

8.79 (7.94-9.42) 

Total alkalinity 

mg CaCOj-U 

4150 (200-54100) 

11795 (300-54100) 

752 (200-2030) 

DOC 

mg-L 

32 (18-82) 

— 

32 (18-82) 

Total P 

Pg'L 1 

1269 (19-15860) 

3863 (243-15860) 

116 (19-273) 

Total dissolved P 

Pg'L J 

762 (14-10866) 

2276 (47-10866) 

89 (14-242) 

Soluble reactive P 

Pg'L' 

454 (3-5816) 

1345 (18-5816) 

57 (3-199) 

Total N 

Pg'L' 

5145 (587-23297) 

12380 (997-23297) 

1930 (587-4379) 

Nitrate + nitrite 

Pg'L' 

31(1-147) 

83 (69-147) 

8 (1-18) 

Dissolved Si 

Pg'L ' 

14 (0.2-93) 

43 (7-93) 

1 ( 0 . 2 - 2 .5) 

Total Fe 

Pg'L ' 

737 (39-8584) 

2308 (141-8584) 

126 (39-574) 

Dissolved Fe 

Pg'L ' 

41 (3-289) 

82 (13-289) 

22 (4-88) 

C:N 


9.6 (6.1-15.1) 

— 

9.6 (6.1-15.1) 

C:P 


231 (79-579) 

— 

231 (79-579) 

N:P 


24 (11-41) 

— 

24 (11-41) 

APA:Chl a 


0.0484 (0.0000-0.4044) 

— 

0.0484 (0.0000-0.4044) 

Conductivity 

mS-crn”' 

9.86 (0.58-72.81) 

18.24 (0.58-72.81) 

6.13 (1.39-16.20) 

Ca 2 * 

mg-L' 

59 (4-338) 

32 (4-127) 

69 (7-338) 

Mg 2 * 

mg-L ! 

179 (0-1729) 

18 (0-65) 

242 (31-1729) 

Na* 

mg-L~' 

1657 (72-17569) 

4123 (72-17569) 

698 (101-2034) 

K* 

mg-L ! 

567 (14-5528) 

1842 (30-5528) 

71 (14-188) 

Cl 

mg-L 

419 (11-4422) 

788 (11—4422) 

255 (18-931) 

so 4 3 - 

mg-L 

2017 (8-10842) 

1186 (8-6138) 

2386 (193-10842) 

Filtered APA rate 

nmol-L MUP-h ' 1 

76.9 (0.0-417.0) 

— 

76.9 (0.0-278.0) 

Unfiltered APA rate 

nmol-L MUP-h ' 1 

31.8 (17.6-82.2) 

— 

31.8 (17.6-82.2) 

a 350 

m' 

23 (2-124) 

69 (25-124) 

8 (2-13) 

a 250 :a 365 


24 (11-71) 

30 (14-71) 

22 (11-45) 

Summer 

Primary production rate 

mg C-m 3 -h ' 

286 (11-1747) 

527 (11-1747) 

125 (15-544) 

Chl a 

Pg'L-' 

118 (3-1613) 

237 (3-1613) 

39 (6-197) 

7 

max 

m 

3.7 (1.0-13.0) 

1.6 ( 1 . 0 - 2 . 6 ) 

5.1 (1.0-13.0) 

Temperature 

°C 

19 (15-24) 

20 (18-24) 

18 (15-23) 

pH 


9.36 (8.41-10.66) 

9.91 (8.88-10.66) 

9.00 (8.41-10.14) 

Total alkalinity 

mg CaCOj-L ' 1 

3725 (190-38000) 

8121 (231-38000) 

794 (190-2550) 

DOC 

mg-L' 

161 (14-1512) 

358 (48-1512) 

29 (14-91) 

Total P 

Pg'L - 1 

4922 (25-71622) 

12087 (89-71622) 

145 (25-329) 

Total dissolved P 

Pg'L - 1 

810 (11-12328) 

1873 (19-12328) 

102 (11-315) 

Soluble reactive P 

Pg'L - 1 

475 (7-5874) 

1085 (7-5874) 

68 (7-271) 

Total N 

Pg'L - 1 

6042 (708-23020) 

11246 (1689-23020) 

2573 (708-3500) 

Nitrate + nitrite 

Pg'L - 1 

8 (0-55) 

11 (0-55) 

5 (1-18) 

Dissolved Si 

Pg'L - 1 

23 (1-156) 

43 (2-156) 

9(1-26) 

Total Fe 

Pg'L - 1 

452 (24-4119) 

925 (24-4119) 

136 (24-516) 

Dissolved Fe 

Pg'L - 1 

31 (0-190) 

26 ( 0 - 112 ) 

34 (3-190) 

C:N 


12.2 (6.3-33.6) 

17.2 (8.5-33.6) 

8.9 (6.3-18.7) 

C:P 


2008 (131-8803) 

4581(432-8803) 

292 (131-1000) 

N:P 


129 (16-705) 

273 (52-705) 

32 (16-53) 

APA:Chl a 


0.0761 (0.00002-0.9667) 

0.1264 (0.00002-0.9667) 

0.0426 (0.00002-0.3309) 

Conductivity 

mS-crn”' 

10.39 (0.53-82.10) 

16.30 (0.53-82.10) 

6.45 (1.51-16.60) 

Ca 2 * 

mg-L' 

45 (4-350) 

11 (4-22.1) 

67 (5-350) 

Mg 2 * 

mg-L' 

168 (0-1889) 

20 (0-62) 

266 (32-1889) 

Na* 

mg-L' 

1753 (74-21073) 

3088 (74-21073) 

862 (101-2716) 

K* 

mg-L' 

743 (15-7954) 

1728 (35-7954) 

86 (15-268) 

Cl 

mg-L' 

473 (13-6394) 

861 (13-6394) 

235 (17-876) 

SO4 2 - 

mg-L' 

2223 (5-12728) 

1459 (5-9176) 

2690 (231-12728) 

Filtered APA rate 

nmol-L ' MUP-h ' 1 

2180.1 (2.5-58734.1) 

5267.0 (2.5-58734.1) 

122.1 (2.8-904.8) 

Unfiltered APA rate 

nmol-L 1 MUP-h ' 1 

4083.1 (5.8-98604.5) 

8979.0 (6.8-98604.5) 

819.1 (5.8-9848.1) 

a 350 

m-' 

20 (2-79) 

38 (11-79) 

8 ( 2 - 22 ) 

a 250 :a 365 


26 (11-87) 

37 (14-87) 

20 (11-38) 


Note: Values are expressed as mean with range in parentheses. Seston ratios are reported on a molar basis. CGP, central Great Plains; NGP, 
northern Great Plains; —, parameter not measured for these sites. APA, alkaline phosphatase activity; MUP, 4-methylumbelliferyl phosphate; 
a 2 so, CI 350 , and a 365 , absoiptivity at 250,350, and 365 nm, respectively. 
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Table 2. R 2 values for simple regressions of rates of primary productivity 
versus selected nutrient parameters. 


Parameter 

Spring 


Summer 


All lakes 

NGP only 

All lakes 

NGP only 

Total P 

0.306" 

0.034 

0.316" 

0.198 

Dissolved P 

0.072 

0.095 

0.011 

0.099 

SRP 

0.030 

0.150 

0.062 

0.055 

Total N 

0.202" 

0.008 

0.168" 

0.238" 

Nitrate + nitrite 

0.561" 

0.103 

0.077 

0.001 


Note: NGP, northern Great Plains; SRP, soluble reactive phosphorus. 
“p < 0.05 


Table 3. AAIC c values describe model fit, with smaller values indicating better fit: (a) summary of primary productivity models devel¬ 
oped with Akaike’s information criterion (AIC) for 2004 survey data; ( b ) summary of top predictor variables from AIC models of 2004 
survey data. 


(a) Summary of primary productivity models developed with AIC for 2004 survey data. 

Model 

Model predictors 

AAICc 

CO 

K 

R 2 

Spring - all lakes (44)* 

log NO 3 (+), C:N (+), log DFE (-), log 0250:0365 (-) 

0.00 

0.086 

4 

0.795 


log NO 3 (+), C:N (+), log DFE (-) 

0.72 

0.060 

3 

0.756 


logCa (-), log SO 4 (-), logTN (+), log DFE (-) 

1.01 

0.052 

4 

0.786 

Spring - NGP lakes only (3)* 

log Cond (-), logCa (-), log TP (+), log DFE (-) 

0.00 

0.652 

4 

0.753 


log Aik (-), logCa (-), log TP (+), log DFE (-) 

2.40 

0.196 

4 

0.717 


logCa (-), log SO 4 (-), log TP (+), log DFE (-) 

2.91 

0.152 

4 

0.709 

Summer - all lakes (58)* 

log Na (-), log Cl (+), C:P (+), log DFE (+) 

0.00 

0.060 

4 

0.715 


Temp (-), log Na (-), log Cl (+), C:P (+) 

0.03 

0.059 

4 

0.714 


log K (-), log TP (+), C:P (+), log DFE (+) 

0.17 

0.055 

4 

0.713 

Summer - NGP lakes only (17)* 

log TP (+), log SRP (-), NP (+), log DFE (+) 

0.00 

0.202 

4 

0.730 


log TP (+), log SRP (-), NP (+), log TFE (+) 

0.65 

0.147 

4 

0.720 


log TP (+), log SRP (-), NP (+) 

1.45 

0.098 

3 

0.636 


( b ) Summary of top predictor variables from AIC models of 2004 survey data. 


95% Cl 

Model-averaged 


Group 

Predictor 

Predictor weights 

regression coefficient 

Lower 

Upper 

Spring - all lakes 

log DFE 

0.94 

-0.794 

-1.239 

-0.349 


log NO 3 

0.81 

0.916 

0.444 

1.387 


C:N 

0.39 

0.071 

0.004 

0.138 


log SO 4 

0.25 

-0.348 

-0.708 

0.012 

Spring - NGP lakes only 

logCa 

1.00 

-1.318 

-2.253 

-0.382 


log DFE 

1.00 

-0.901 

-1.283 

-0.520 


log TP 

1.00 

1.699 

0.916 

2.482 


log Cond 

0.65 

-1.124 

-1.628 

-0.620 

Summer - all lakes 

C:P 

1.00 

0.00015 

0.00009 

0.00021 


log Cl 

0.57 

0.458 

0.151 

0.766 


log K 

0.53 

-0.485 

-0.792 

-0.178 


logNa 

0.47 

-0.604 

-0.948 

-0.261 

Summer - NGP lakes only 

log TP 

0.69 

2.440 

1.216 

3.665 


N:P 

0.68 

0.018 

0.005 

0.031 


log SRP 

0.62 

-1.338 

-2.201 

-0.476 


Note: Akaike weights, denoted as to, describe the probability that the model is the best model generated, and K is the number of estimable parameters 
in the model. Predictor weights indicate the predictive strength of the variable, with values near 1 having the strongest predictive weight. Asterisk (*) 
indicates the number of models with AAIC c < 4, which were used to compute the Akaike weights. Cl, confidence interval; DFE, dissolved Fe; TN, total 
N; Cond, conductivity; TP. total P; Aik, alkalinity; Temp, temperature; SRP, soluble reactive phosphorus; TFE, total Fe. 


only responded to the N treatment (F = 81.3, df = 4,10, p < 
0.001), with no additional increase in the combination treat¬ 
ment of N + P. It should be noted, however, that Free Peo¬ 
ple Lake was classified as a SPR-NP lake based on 2004 


data, but in 2006, the APA: Chi a ratio was below 0.005, 
moving it to the SPR-N group (Table 5). For the SPR-N 
group, predicted to respond to N alone, the rate of produc¬ 
tion only increased in East Devil’s Lake in the N addition 
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Fig. 2. Lake classification for spring 2004 survey data. Nitrate, shown at the top of the figure, was identified as the main splitting variable 
based on an initial regression tree with primary production as the response metric. On the lower left portion of the figure, NGP lakes were 
further examined for P limitation with APA: Chi a ratios identified as the splitting variable based on a classification tree with P-limitation 
category (no P limitation or severe limitation) as the response variable. Simple regressions were performed within each subset of lakes, 
although no regressions were statistically significant (a = 0.05). Abbreviations: CGP, central Great Plains; NGP, northern Great Plains; 
SPR-N, potentially N-limited lakes; SPR-NP, potentially N and P co-limited lakes. 



(F = 96.5, df = 4,10, p < 0.001). Although there were signif¬ 
icant differences across all treatments for Stink Lake (F = 
3.54, df = 4,10, p = 0.048), differences were not large 
enough for Tukey’s post-hoc test to identify homogeneous 
subsets. There appeared to be a bottle effect for this lake, 
with a substantial increase in production rates in all treat¬ 
ments compared with the initial values (for example, 
669 mg C-m Mi 1 in the control compared with the initial 
value of 31.5 mg C-m -3 -h _1 ). 

Summer experimental results 

For the summer experiment, three lakes were tested from 
each of the SUM-N (Alkaline Lake, East Devil’s Lake, and 
Stink Lake) and SUM-P (Clear Lake, Coldwater Lake, and 
George Lake) groups. Rates of primary production differed 
among nutrient treatments in five of the six lakes (Fig. 5; 
Table 4). Of the SUM-N lakes predicted to respond to N, 
production rates in East Devil’s Lake increased in response 
to the N addition (F = 25.4, df = 4,10, p < 0.001). In Alka¬ 
line Lake, rates increased in both the N treatment and the P 


treatment (F = 80.9, df = 4,10, p < 0.001), with even higher 
production rates in the N + P combination. Interannual dif¬ 
ferences may have also been important in Alkaline Lake, as 
this lake was originally classified in the SUM-N group 
based on 2004 data (Table 5). In the summer of 2006, how¬ 
ever, the SRP concentration was undetectable (0.010 pg-L ! 
detection limit), changing its classification to SUM-P. There 
were no significant differences in production rates for Stink 
Lake (F = 1.329, df = 4,10, p = 0.325) due to apparent bot¬ 
tle effects for this lake, as in the spring experiments. Control 
treatment production rates were elevated compared with the 
initial samples (167 mg C-m Ur 1 in the control versus an 
initial rate of 23.1 mg C-m'Lfr 1 ). In the SUM-P lakes, pre¬ 
dicted to be P-limited, Clear Lake had the highest produc¬ 
tion rates in the P treatment (F = 26.3, df = 4,10, p < 
0.001), with no additional response to the N + P treatment. 
Production rates in Coldwater (F = 195.9, df = 4,10, p < 
0.001) and George (F = 92.0, df = 4,10, p < 0.001) lakes 
both increased in the N addition, along with a significant in¬ 
crease in the N + P addition. 
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Fig. 3. Lake classification for summer 2004 survey data. C:P ratios, at the top of the figure, were identified as the main splitting variable 
based on an initial regression tree with primary production as the response metric. On the upper right, high C:P lakes from the CGP region 
showed a nonsignificant relationship between production and NO 3 - (p = 0.073, y = 0.281x + 2.636, df = 1,5, F = 5.12). On the lower left of 
the figure, lakes from the NGP region were further examined for P limitation with a classification tree with APA: Chi a class (no P limita¬ 
tion or severe limitation) as the response metric. Simple regressions were performed within each subset of lakes. The regressions were sta¬ 
tistically significant for the SUM-N (p = 0.043, y = 0.0246.V + 1.219, df = 1,9, F = 5.55) and SUM-P (p = 0.009, y = 1.063* - 0.194, df = 
1,5, F = 17.5) groups. Abbreviations: CGP, central Great Plains; NGP, northern Great Plains; SUM-N, potentially N-limited lakes; SUM-P, 
potentially P-limited lakes. 



Discussion 

Although previous studies of algal production in saline 
lakes have been able to identify limitation patterns within a 
few systems (Bierhuizen and Prepas 1985; Campbell and 
Prepas 1986; Waiser and Robarts 1995), our study investi¬ 
gated drivers of primary production by phytoplankton across 
a diverse set of saline lakes. The combination of predictive 
modeling and experimental tests allowed us to begin to de¬ 
lineate subsets of lakes and determine limitation patterns 
within them. The models and experiments presented here in¬ 
dicate that despite the high total nitrogen and phosphorus 
concentrations in these lakes, phytoplankton in these lakes 
are potentially nutrient-limited. Though many freshwater 
lakes are often phosphorus-limited (Dillon and Rigler 1974; 
Wetzel 2001), our results reveal widespread nitrogen limita¬ 


tion across these prairie saline lakes. A subset of these lakes 
is also secondarily limited by P, as indicated by high 
APA: Chi a ratios and large additional increases in algal 
production in the experimental N + P treatments. Although 
there is evidence for nitrogen limitation in single saline sys¬ 
tems (Reuter et al. 1993), our study provides information on 
nutrient limitation patterns across a large region in central 
North America with a high density of subsaline and saline 
lakes. These data support observations that nitrogen may 
play a larger and more widespread role in nutrient limitation 
patterns in lake ecosystems than previously thought (Smith 
1982; Elser et al. 1990; Lewis and Wurtsbaugh 2008). 

Multiple statistical techniques enabled us to begin to dis¬ 
cern these productivity patterns in prairie saline lakes, and 
results from these various methods were generally in agree- 
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Fig. 4. The results from a nutrient enrichment bioassay for spring rates of primary production that experimentally tested the lake classifica¬ 
tion groups: ( a ) SPR-NP and ( b ) SPR-N. Abbreviations: Aik, Alkaline Lake; CW, Coldwater Lake; FP, Free People Lake; ED, East Devil’s 
Lake; Stk, Stink Lake; Prim prod, rate of primary production. 



Treatment 


Table 4. Summary of Tukey’s post-hoc analyses for spring and summer 2006 experiments. 


Treatment 

Spring 




Summer 





Aik* 

CW* 

FP* 

ED* Stkt 

Aik* 

ED* 

Stk* 

Clear* 

CW* 

Grg* 

C 

1 

1 

1 

1 — 

1 

1 

— 

1 

1 

1 

Fe 

1 

1 

1 

1 — 

1, 2 

1 

— 

1 

1 

1 

P 

1, 2 

1 

1 

1 — 

2 

1 

— 

2 

1 

1 

N 

2 

2 

2 

2 — 

2 

2 

— 

1 

2 

2 

N + P 

3 

3 

2 

2 — 

3 

2 

— 

2 

3 

3 


Note: Numbers in columns indicate treatments within each lake with responses that are significantly different from each other (a = 0.05). Nutrient treat¬ 
ments: C, control; Fe, Fe addition; P, P addition; N, N addition; N + P, combined N and P additions. Lake abbreviations: Aik, Alkaline Lake; CW, Cold- 
water Lake; FP, Free Peoples Lake; ED, East Devil’s Lake; Stk, Stink Lake; Clear, Clear Lake; Grg, George Lake. 

*p < 0.001 for analysis of variance. 

1 p = 0.044 for analysis of variance, but post-hoc analysis did not detect different subsets of lakes, 
y? = 0.325 for analysis of variance. Post-hoc analysis was not performed on these data. 


Table 5. Summary of changes in clas¬ 
sification groups from 2004 to 2006. 



Survey class 

Lake 

2004 

2006 

Spring 

East Devil’s 

SPR-N 

SPR-N 

Stink 

SPR-N 

SPR-N 

Free Peoples 

SPR-NP 

SPR-N 

Alkaline 

SPR-NP 

SPR-NP 

Coldwater 

SPR-NP 

SPR-NP 

Summer 

East Devil’s 

SUM-N 

SUM-N 

Stink 

SUM-N 

SUM-N 

Alkaline 

SUM-N 

SUM-P 

Coldwater 

SUM-P 

SUM-P 

George 

SUM-P 

SUM-P 

Clear 

SUM-P 

SUM-P 


Note: Lakes given in italic type indicate 
changed from predicted 2004 classification 
based on 2006 data. 


ment. Simple and multiple regressions identified some po¬ 
tential limiting factors (such as NO 3 - in the spring for all 
lakes), but did not suggest singular seasonal controls across 
the entire suite of lakes. In addition, with covariation among 
some key parameters in these systems (e.g., correlations be¬ 
tween logTN and log TP were R 2 = 0.88 in spring and R 2 = 
0.64 in summer), it was not possible to identify whether one 
or both nutrients were actually limiting. Therefore, using ad¬ 
ditional statistical techniques, we developed more complex 
models to decipher drivers of productivity patterns across 
these diverse systems. AIC revealed a complex multivariate 
story of predictors of primary production, identifying several 
different predictive models and varying parameter weights. 
The CART analyses identified a similar group of variables 
to those from regressions and AIC but produced a simplified 
hierarchical ordering of these variables and thus served as a 
useful exploratory tool for examining broad trends across an 
ecologically complicated data set, as suggested by De’ath 
and Fabricius (2000). We recognize that the sample size for 
this study was relatively small, particularly after separating 
the lakes into subsets, and that more lakes would strengthen 
these models. However, even with a small sample size, the 
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Fig. 5. The results from a nutrient enrichment bioassay for summer rates of primary production that experimentally tested the lake classifi¬ 
cation groups: (a) SUM-N and ( b ) SUM-P. Abbreviations: Aik, Alkaline Lake; ED, East Devil’s Lake; Stk, Stink Lake; Clear, Clear Lake; 
CW, Coldwater Lake; Grg, George Lake; Prim prod, rate of primary production. 



C Fe N P NxP C Fe N P NxP 

Treatment 


CART models helped to clarify trends within lake groups. 
Classifications developed from key variables identified with 
these techniques were easily tested via experiments in 2006. 
In the spring, the CART models correctly predicted N limi¬ 
tation in the SPR-N and SPR-NP lakes and secondary P lim¬ 
itation in the SPR-NP lakes. In the summer, lakes in the 
SUM-N group responded to N as predicted, except for one 
lake that changed classification groups in 2006 and was co¬ 
limited by N and P. We predicted an increase in production 
in P treatments in the SUM-P lakes, and although there was 
a response to P, two of the three lakes were mainly limited 
by N and secondarily by P. Despite this misidentification 
and generally small sample size, the CART models as a 
whole enabled us to identify limitation patterns across prai¬ 
rie saline lakes more clearly than in prior studies. 

Although ambient nutrient concentrations were high, ex¬ 
perimental nutrient additions stimulated primary production 
by phytoplankton. One explanation for this is the reduced 
bioavailability of nutrients in these lakes due to complexa- 
tion with DOM (Hessen and Tranvik 1998; Findlay and Sin- 
sabaugh 2003). As direct measurements of nutrient 
concentrations did not predict rates of primary production, 
physiological indicators of nutrient limitation were very im¬ 
portant to the success of the lake classification system, with 
APA: Chi a serving as a reliable tool in discerning patterns 
of phosphorus limitation among the lakes. Unfortunately, a 
similarly reliable indicator of nitrogen limitation in these 
lakes is not apparent. Although N 2 -fixation rates and hetero¬ 
cyst formation can indicate nitrogen limitation for cyanobac¬ 
teria, few heterocysts were found in phytoplankton samples 
from the 2004 sampling (personal observation). C:N ratios 
of particulate material can be indicative of nitrogen limita¬ 
tion in other aquatic systems, with low C:N ratios suggest¬ 
ing that algae are nitrogen replete (Wetzel 2001); however, 
C:N ratios never emerged as significant in our analyses. Fur¬ 
thermore, in many lakes that responded positively to N addi¬ 
tions, initial seston C:N values were low (many <9), which 
clearly did not reflect the physiological status of the algae. 


This observation, in conjunction with the extremely high 
seston C:P ratios in these lakes, raises the question of what 
these C ratios really reflect in these high DOM lakes. The 
development of additional physiological metrics would 
greatly improve our ability to assess N limitation patterns in 
these prairie lakes. 

Our results indicate broad nitrogen limitation in lakes 
across this semi-arid region; other studies have also found 
N limitation in single saline lakes (Cloern et al. 1983; 
Herbst and Bradley 1989; Reuter et al. 1993). Terrestrial 
sources of nitrogen may be less abundant in grasslands than 
in forested areas, as humid grasslands tend to be nitrogen- 
limited with low soil N availability (Wedin 1995). These 
systems may be further driven to N limitation by various 
feedbacks, as dry conditions in grasslands slow N minerali¬ 
zation rates of already low quality (high C:N) litter (Wedin 
1995). Other important feedbacks to consider in these prairie 
ecosystems include grazing and fire effects, as grazing tends 
to recycle N (Holland et al.1992) and fire tends to volatilize 
N (Ojima et al. 1994). In contrast to wetter temperate re¬ 
gions, reduced precipitation in grasslands can also decrease 
nutrient inputs from terrestrial runoff. 

Nitrogen limitation in prairie saline lakes may also be 
caused by different physical and chemical conditions in 
these systems, which lead to N loss from the system and P 
release from sediments. Physically, lakes in the CGP and 
NGP tend to be shallow and turbulent, and the resulting var¬ 
iable mixing patterns could lead to alternating aerobic and 
anaerobic conditions at the lake bottom. These varying con¬ 
ditions allow for high rates of denitrification, a common 
pathway for nitrogen loss in shallow systems (Jensen and 
Andersen 1990). Though aerobic conditions typically cause 
P to precipitate out of the water column due to binding with 
Fe, turbulence also resuspends sediment and may cause 
overall P release if conditions vary between aerobic and 
anaerobic (Spndergaard et al. 1992). Chemically, Fe-P bind¬ 
ing is reduced at high pH (Lijklema 1977), and all study 
lakes had pH greater than 8.3. There is also evidence for al- 
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tered P cycling under high sulfate concentrations, as the sul¬ 
fate tends to bind Fe and thereby reduce Fe-P binding 
(Caraco et al. 1989; Curtis 1989). Saline lakes, particularly 
those in the NGP, can have high sulfate concentrations, and 
although ion concentrations did not directly correlate with 
primary production rates, sulfate may indirectly influence 
nutrient availability in these systems. 

Although Fe was a top predictor variable for primary pro¬ 
duction in AIC and multiple regression models, the CART 
analyses did not identify Fe as important, and there was lit¬ 
tle evidence for Fe limitation in experimental results. In the 
spring experiments, a Chroococcus sp. in Stink Lake did in¬ 
crease substantially in biovolume in the Fe treatment (per¬ 
sonal observation), suggesting that Fe was biologically 
available for uptake in this lake and season. Our results 
show that limitation by two or possibly more nutrients could 
be very important in controlling production in prairie lakes, 
and previous studies have indicated the potential of Fe and 
other trace metals (Marino et al. 1990; Evans and Prepas 
1997) to act as limiting factors in prairie saline lakes. Future 
studies should consider the effects of not only single nutrient 
additions, but also those of combination treatments as well. 

Lakes in the CGP were generally more productive than 
those of the NGP and, in the spring, had higher nitrate con¬ 
centrations. Precipitation gradients are similar in both re¬ 
gions, but the higher N and hence productivity of the CGP 
lakes may be due to the higher temperatures in this region, 
which could increase terrestrial N mineralization rates, as 
well as N 2 fixation rates by cyanobacteria (Macarelli and 
Wurtsbaugh 2006; Scott et al. 2008). Land-use patterns also 
differ across these two regions, with the CGP serving as cat¬ 
tle rangeland and the NGP largely serving as arable agricul¬ 
tural land, some of which has been abandoned during recent 
decades. These differences in land use, as well as local geol¬ 
ogy as described above, may contribute to the higher pro¬ 
ductivity levels in the CGP lakes. 

Shallow lakes and wetlands in prairie regions provide im¬ 
portant ecosystem services, including carbon sequestration, 
wildlife habitat, and recreation (Batt et al. 1989; Euliss et 
al. 2006). These ecosystem services all involve algal produc¬ 
tion and standing crops, and hence a greater understanding 
of the factors that control primary production across these 
lakes is a critical aspect for ensuring the sustainability of 
these services. Our analyses suggest that patterns of primary 
production in prairie saline lakes of the central and northern 
Great Plains were complex. Although no trends were consis¬ 
tent across all saline lakes, likely due to geochemical and 
physical differences among these varied systems, explora¬ 
tory classification models generally were successful in iden¬ 
tifying productivity patterns within lake subsets. Nitrogen 
limitation was widespread in these systems, especially in 
spring, along with secondary limitation by phosphorus in 
some cases. The CART exploratory classification technique 
shows promise for determining primary production patterns 
across a complicated data set, and the resultant classification 
could be refined with a larger data set as well as with better 
physiological tools to assess limitation by nitrogen and pos¬ 
sibly other nutrients. Prairie saline lakes are globally impor¬ 
tant systems, and thus, understanding the ecology of these 
complex systems is an essential component of a more com¬ 


prehensive understanding of lake ecosystem function and 
structure. 
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